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Abstract. A discretization of the peakons lattice is introduced, belonging to 
the same hierarchy as the continuous-time system. The construction examplifies 
the general scheme for integrable discretization of systems on Lie algebras with 
r-matrix Poisson brackets. An initial value problem for the difference equations is 
solved in terms of a factorization problem in a group. Interpolating Hamiltonian 
flow is found. A variational (Lagrangian) formulation is also given. 



1 Introduction 



The subject of integrable symplectic maps received in the recent years a 
considerable attention. Given an integrable system of ordinary differential 
equations with such attributes as Lax pair, r-matrix and so on, one would like 
to construct its difference approximation, desirably also with a (discrete-time 
analog of) Lax pair, r-matrix etc. Recent years brought us several successful 
examples of such a construction [1-9]. 

Recently there appeared for the first time examples , where the Lax 
matrix of the discrete-time approximation coincides with the Lax matrix of 
the continuous-time system, so that the discrete-time system belongs to the 
same integrable hierarchy as the underlying continuous-time one (systems of 
Calogero-Moser type). This led us to the formulation of a general recipe for 
producing discretizations of this type, which was applied to the two favorite 
examples of finite dimensional integrable systems of the classical mechanics 
- the Toda lattice and the relativistic Toda lattice 0,0. 

In the present Letter we want to describe an application of this scheme to 
another remarkable system, discovered recently by Camassa and Holm [pi]|]. 
Such an application was made possible thanks to an r-matrix interpretation 
given to this system by Ragnisco and Bruschi |TT . 



2 Continuous— time peakons lattice 

The peakons lattice has appeared in |TD| as a system describing positions and 
velocities of "peaked solitons" for a certain integrable shallow water equation. 
The Hamiltonian function of the peakons lattice reads: 

1 ^ 

Hix,p) = PjPkexp{-\xk-Xj\). 

It can be proved |10[ that, generally speaking, the peakons retain their or- 
dering, so that one can assume x\ > X2 > ... > xn, and the Hamiltonian 
function takes the form 



TV 



H{x,p) = -J2pI+ Y1 PjPk exp(xfc 

A,-l l<j<k<N 



(2.1) 
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The canonical equations of motion read: 

^k = J2 Pj exp(a;fc - Xj) +Pk + Yl Pj exp(xj - Xk), (2.2) 

j<k j>k 

— = - J2Pj exp(xfc - Xj) + J2 Pj exp(xj - Xk). (2.3) 

P!<^ j<k j>k 



One can consider (|2.2| ) as a linear system for the impulses pk- Remarkably 
enough, the matrix of this linear system has a tridiagonal inverse, so that pk 
depends only on Xj, Xj with j = k — l.k^k + 1: 

1 -exp(2(xfc+i -a;fc_i)) 

Pk = Xk] 



[1 - exp(2(xfc+i - Xk))][l - exp{2{xk - Xk-i))] 

exp{xk+i-Xk) . exp(xfe-Xfe_i) 

- Xk+i- 77-^ - Xk-i- 7-7 TT-. [^A) 

1 - exp(2(xfc+i - Xk)) 1 - exp(2(xfe - Xk-i)) 

For k = 1 and k = N little modifications are necessary, which can be taken 
into account by imposing the boundary conditions 

Xq = 00, xn+1 = —00. (2.5) 

The formula ( |2.4| ) implies that the Lagrangian function for the peakons lattice 
is local (with only nearest-neighbour contributions): 

x) = -J2 xl- ^ ~ exp(2(x,.+i - Xk-i)) 



2 ^1 [1 - exp(2(xfc+i - Xk))][l - exp(2(xfc - Xk-i))] 



N-l 



^ . . exp(Xfc+i - Xkj 1^ „s 

XkXk+X- 77^7 TT, (2.6) 

f-: 1 - exn 2 Uct_n - 



k=\ 



exp(2(a;fc+i - x^))' 



The Lax pair formulation with the matrices from gl{N) for the system 



(|2.2|) - (p.3|) was also given in 

t = [T, A] (2.7) 

with the Lax matrix 

^ /I \ 

T{x,p) = Y^PkEkk + VPkPj exp I -{xk - Xj) ] {Ekj + Ejk), (2.8) 
k=i j<k ^ 
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and ^ ^ 

Mx,p) = 7^J2VPkPj^^P [T^ixk-Xj)] {Ekj-Ejk)- (2.9) 

An r-matrix interpretation of this result was given in [0. To formulate 
it, let (■, ■) denote the standard scalar product in gl{N): {X,Y) = tr(Xy), 
and let the gradient Vip(T) of a smooth function ip on gl{N) be defined by 



VX G gl{N). 



e=0 

Let R{X) be a linear operator on gl{N), defined as 

strictly lower triangular part (X) — strictly upper triangular part (X), 

This operator satisfies the classical modified Yang-Baxter equation, so that 
the following Poisson bracket on gl{N) is defined: 

W^i^UT) = i(T, [i?(v^(r)), v^(T)] + [V¥.(T),i?(v^(r))]). 

Now one of the main results of 0] states that the set of the matrices T{x,p) 
from forms a Poisson suhmanifold in gl{N) equipped with the above 

Poisson bracket. Hamiltonian equations generated in this bracket by a con- 
jugation invariant Hamiltonian function (p{T) read: 



T 

where 



T, ^i?(Vy.(T))] = [t, 7r±(V^(r))] , (2.10) 



7r+ = -(i? + /), 7r_ = -(i? — /), so that 7r+ — 7r_ = /. 

This explains the formulas ( p.7|) -( pl9D , since the Hamiltonian H{x,p) of the 
peakons lattice is nothing else then ip(T) = |tr(T^), so that VfiT) = T. 



3 Discrete— time peakons lattice 

A general recipe for obtaining integrable time discretizations for equations of 
the form ( |2.1CI| ) was given in 0,0. Discrete time systems obtained by this 
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recipe share with the underlying continuous-time ones the Lax matrix, and 
hence the integrals of motion and the integrability property. The formula 
reads: 

T{t + h) = U^Hl + hVv{T{t)))T{t)uJl + hV^{T{t))) (3.1) 



Here the factorization Y = Il^{Y)Ilz^{Y) in the group GL{N) corresponds 
canonically to the additive decomposition X = 7r+(X) — 7r_(X) in its Lie 
algebra gl{N), and in our case is characterized by the following conditions: 
n_|.(F) (n_(F)) is a nondegenerate lower triangular (resp. upper triangular) 
matrix, and the diagonals of these matrices are mutually inverse. 

The aim of the present Letter is to apply this recipe to the peakons lattice, 
i.e. to put in the equation ( |3.1| ) the Lax matrix ( p.8| ) and Vip{T) = T. 

Theorem 1. The equation 

T{t + h) = n;i(/ + /iT(t))T(t)n+(/ + /iT(t)) (3.2) 

serves as a Lax representation for the following map Q in the space R'^'^{x,p}: 
exp{xk - Xk) = + hpk + /i^pjexp(xj - Xk), (3.3) 

!^ = + 1 J^hl R E P3 exp(a^i - Xk) ■ (3.4) 
Pk J- + iiPkPk jyk 

Here the quantities (3k are defined recurrently by the relation 

n . I X , eXp(Xfc+l - Xk) , . 

Pfe+i = 1 - exp(a;fc+i - Xfc) H j , (3.5) 

hpk + -r 

Pk 

with an initial condition 

A = i- 

Remark 1. A simple induction shows that the finite continued fractions 
Pk have by /i — the following asymptotics: 

= 1 - /i^pj exp(xfe -Xj) + 0{h^). (3.6) 

j<k 



^We use the tilde to denote the time h shift, so that Xk = Xk{t + h), if Xk = Xk{t). 
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Hence the discrete time equations ( p.3| ), ( p.4| ) in fact approximate the con- 
tinuous time ones (|2.2|) , (|2.3|) . 

Remark 2. Due to the general theory (cf. ||^, 0), the initial value 
problem for the discrete time peakons lattice can be solved in terms of 
the Cholesky factorization of the matrix (J + /iT(0))", and the interpolat- 
ing Hamiltonian for this system is given by tr($(T)), where 

= t log(l + hr])dr] = + 0{h). 
Jo 2 

Proof of the Theorem 1. Due to the symmetry of the matrix T, the 
n_|_ni^ factorization for / + hT takes the form of the so-called Cholesky 
factorization: 

I + hT = LL^, (3.7) 

where 

L = Ii+{I + hT) = Y,h^Euj (3.8) 

k>j 

is the (uniquely defined) lower triangular matrix with positive diagonal en- 
tries. With the help of (p.7|), (p.8|) the equation (p.2|) can be rewritten as 



I + hT = L^L. (3.9) 



It turns out that the special structure of the matrix T p.^ ) allows for an 
almost explicit computation of the factor L. 

Lemma. 

i) For k > j 



/fe+i fc+i = 1 + hpk+i - exp(xfc+i - Xfc)^^ ( 7^ - 1 + ^Pk] ■ (3.11) 

The proof of this lemma is relegated to the end of this section. The 
Theorem 1 simply follows from this lemma. Indeed, setting in ( |3.11| ) 



ll = l + hpkPk, (3.12) 
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we immediately arrive at the recurrent equation ( |3.5| ) for Pk- The following 
relation, equivalent to ( p.l2|) , appears to be also useful: 



hk — - — 



IkkJ 1 + hpkPk 



(3.13) 



Substituting now (p.8|) into ( p.9|) , we get: 



l + hpk = Y^ lik, h^JpkPj exp (^(xfe - Xj)) = likhj for k > j. (3.14) 



i>k 



i>k 



With the help of ( |3.10D the first equation in ( ^.14| ) may be re-written as 



hpk = Ilk - + (hk - J-) exp{xi - Xk). 



i>k 



Pk 



(3.15) 



Now ( p.4|) follows immediately, upon application of ( |3.12| ), (|3.13| ). 

The right-hand side of the second equation in (|3.14|) may be brought with 
the help of ( |3.10[ ) into the form 



fPk /I, 
— exp (-(xfc 



\ ''jj / 



kk+ [hk - Y-)J2~ exp(xi - Xk) 

l-kk/ iyk Pk 



Because of ( p.l5|) this is equal to 

[Pk /I. 
'—exp {^-{xk -. 



Pj 



hk — - — 



kk 



Hence we arrive at the conclusion that the quantity 



exp(Xfc - Xk) 
PkPk 



kk 



kk 



does not depend on k. Setting it equal to /i^, we get upon use of ( |3.13D : 

Pk 1 + hpk(3k 



expi^Xfc - Xk) 



Pk PI 



(3.16) 
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This, together with ( p.4| ) is equivalent to ( |3.3| ). The Theorem 1 is proved. 

Proof of the lemma. The statement i) may be proved by induction in 
j. Suppose that it holds for l^i with k > i for all i < j. From the definition 
of the matrix L ( |3.7| ) we have: 

J2 = 1 + hpf, Ikilji = hy/pkPj exp (^{xk - Xj)\ for k > j. (3.17) 



Using the base of induction, i.e. the expressions ( p.lOj ) for l^i, Iji with i < j, 
we deduce from the first equation in ( p.l7| ): 



l+Pj 



Kj 



exp(Xj — X 
Pi 



I'i'i 



(3.18) 



and from the second one: 



Ikjljj = ^/WPj exp {-{xk-Xj)] h-J2 



exp(a;j — x 



i<j 



Pi 



Comparing the right-hand side of the last equation with ( p.l8|) , we obtain: 



Pj 



This is identical with (|3.1CI|) , which is hence proved for all k > j. 
The statement ii) follows now from the formula ( |3.18| ): 

,2 1 A _ exp(a;fc+i - Xi) f 1 ^ ^ 
'fc+i,fc+i - J- - n-Pk+i - -Pk+i 2^ I ~ T 

i<k Pi ^ 



-pfe+iexp(xA,.+i -Xk) 



Pk V ffcfc/ i<fc 



exp^Xfc - Xj 
Pi 



[using again ( |3.18| ) to express the last sum) 



Pk+i 
Pk 



exp(a;fc+i - Xk) 



1 \2 

kk-^] - {Ilk - 1 - hpk) 

f-kky 



and ( p. Ill ) follows. The lemma is proved. 
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4 Lagrangian formulation 

Recall that in the continuous time case, considering the evolution equation 
(|2.2|) as a linear system for pk-, we got a "tridiagonal" expression (p.4|) . In 
the discrete time case, the corresponding evolution equation (|3.3|) is highly 
nonlinear, due to the continued fraction expression ( p.5|) for (3k- Nevertheless, 
it turns out to be possible to obtain "tridiagonal" expressions for the impulses 
Pk in the dicrete time case, also. 

In order to formulate the corresponding result, it will be convenient to 
introduce following short-hand notations: 

1 1 
Afc = -(exp(-Xfc+i) -exp(-Xfc)), 5^ = -(exp(xfc) - exp(xfc+i)). (4.1) 

Theorem 2. There hold following relations: 



hPk = ^(exp(a;fc-i 



Xk) - exp(Xfe+i - Xk 



+ exp(-Xfc)SfeJl + -^— -exp(-Xfc)5fe_.iWl + -- — ^- — ; (4.2) 

V ^fc-Dfc V ^fc-l-Dfc-l 



hPk = ^{exp{xk - Xk+i) - exp(xfc - Xk-\, 



- exp(xfc)ylfcWl + + exp(xfc)Afc_i Jl + - — — . (4.3) 

V ^fc-Dfc V ^fe-l-Dfc-l 

This Theorem will be proved at the end of this section. Several remarks 
are now in turn. 

Remark 3. For k = \ and k = N minor correction in the formulas 



([4.3| ) are necessary, due to the boundary conditions (p.5|). The corresponding 



formulas read: 



hpi = -1 + ^(^exp(xi - Xi) - exp(x2 - Xi)) + exp(-a;i)5i Jl + — ^ 



hpi = ^(^exp(a;i -xa) +exp(xi -xi)) -exp(a;i)AiWl + — i 
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hpN = ^(^exp{xN~i-XN)+exp{xN-XN)) -exp{-XN)BN-iJl + 



An-iBn-i 



hpN = -l + ^(^exp{xN -xn) -exp{xN -xn-i)) +exp{xN)AN_iJl + 



N-l 



Remark 4. The expressions (|4.2| ), (4.3) serve as finite-difference approx- 
imations to (P^). However, as opposed to the analogous statements in [|], 
[P], this is not immediately obvious, but rather requires for some (straight- 
forward) calculations. 

Now we are in a position to give a Lagrangian formulation of the discrete 
time peakons lattice, i.e. to represent it in the form 

d(K{x{t + h),x{t)) + K{x{t),x{t-h)))/dxk{t) =0. (4.4) 

The key point is the representation of the momenta in terms of the La- 
grangian function [|I^: 

Pk = -dA{x,x)/dxk, Pk = dA{x,x)/dxk. (4.5) 



Identifying these expressions with ( [4.2|) , (|4.3|) , we see that the following state- 
ment holds. 

Theorem 3. The discrete time peakons lattice is a Lagrangian system 
with the Lagrangian function 

A{x, x) = Xi + - exp(5;i — xi) + - exp{xN — xn) — xn 



+ 2 51 (exp(5fc - Xfc+i) - exp(xfc+i - Xfc)) - ^ ^(Afc5fc), (4.6) 



fc=i fc=i 

where 



vl>(e) = 2 /\/l + irfr^. 
JO V ri 



r] 

Remark 5. The function ( [4.6| ) serves as a rather non-trivial approxima- 
tion to ( p.6|) . This can be seen after some straightforward, though tediuos 
calculations. 



9 



Proof of the Theorem 2. Subtract from the equation ( p.3| ) the analo- 
gous equation for the subscript k + 1 muhiphed by exp{xk+i — Xk): 



exp(xfc+i - Xk) = exp(xfe - x^) - exp(xfc+i - Xk). 



Using the recurrence ( ^.51 ) to express Pk+i, we end up after some simple 
calculations with the following quadratic equation: 

-5- + hpk] -2eyip{-Xk)Bk ( -3- + hpA - exp(-2a;fc)^ = 0. 



Its (positive) solution: 



In order to exclude f3k from this expression, we use ( ^.51 ) in the form 

/ \ 

exp(-a;fe_i) 



(4.7) 



(3k = exp{xk) 



2Ak-i + 



V 



- — + hpk-1 I 
Pk-i I 



Using in the right-hand side ( [4.7|) , we get 

^fc = exp(xfc)Afc_i Wl + 



1 



Ak-\Bk-\ 



or 



^ = exp(-Xfc)Bfe_i Ml 



(4.8) 
(4.9) 



Ak-iBk-i J 

Now (|]2D follows directly from (^ ) and ( ^ . 

Finally, to prove (|4.3| ) we use ( p.l6|) , re-written in the form 

h~ h f~ I \ (3kexp{-Xk) 
hpk = hpkexp{xk + Xk)-rj • 

[ — + /ipfc j exp{xk) 

Putting (^4.2| ), ( [4.7|) , (|4.^ ) into this formula, we arrive after some manipula- 
tions at ([4.3|). The Theorem 2 is proved. 
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5 Conclusion 



A new application of a general scheme for producing integrable discretizations 
for integrable Hamiltonian flows is described in the present Letter. Advan- 
tages of this approach are rather obviuos: it is, in principle, applicable in a 
standartized way to every system admitting an r-matrix formulation, at least 
with a constant r-matrix satisfying the modified Yang-Baxter equation. We 
shall demonstrate elsewhere that the discrete time systems from 0, with 
dynamical r-matrices may be also included into this framework. We hope 
also to report on numerous other applications of this approach in the future. 

The drawback of this scheme is also obvious to any expert in this field. 
Namely, some of the most beautiful discretizations do not live on the same 
r-matrix orbits as their continuous time counterparts [Q, [|], @], |^, and 
there seems to exist no way of a priori identifying the correct orbit for nice 
discretizations. However, we hope that continuing to collect examples will 
someday bring some light to this intriguing problem. 

The research of the author is financially supported by the DFG (Deutsche 
Forschungsgemeinshaft) . 
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